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^ ' Abstract 

o , 

S^ ' Abstract: We study the infrared behaviour of the two-dimensional Euchdean 0{N) 

CN ■ nonhnear a-Model with a suitable ultraviolet cutoff. It is proven that for a sufficiently 

Q>^ ' large (but finite!) number N of field components the model is massive and thus has 

^^ ' exponentially decaying correlation functions. We use a representation of the model 

with an interpolating bosonic field. This permits to analyse the infrared behaviour 
without any intermediate breaking of 0(A)-symmetry. The proof is simpler than that 
of the corresponding result for the Gross-Neveu- Model [1]. 



X ■ 1 Introduction 



We want to study the infrared behaviour of the two-dimensional Euclidean nonlinear cr-model 
[2] which is formally given in terms of the Lagrangian 

c = ^{{d<pr + ^{<p^-m. (1) 

Here the constant K is assumed to be of order 1, whereas we assume N » 1 , for A see 
below. is a real-valued A^- (flavour-) component bosonic field in the fundamental (vector) 
representation of 0{N). The minimum of C is thus situated at 0^ = 1, where the value 1 
may be changed by rescaling the field variable. The ultraviolet (UV) cutoff as well as more 
precise statements on the lower bound for A^ will be specified later. As regards A, its value 
should not be much larger than 1, because otherwise the generated mass m approaches the 
UV cutoff, see below (20). If it is much smaller than 1, on the other hand, the effective 
energy range of the UV cutoff model becomes large and therefore the bounds, which involve 



factors of exp(47r/A), deteriorate. The convergence proof then requires larger values of A^. 
In the full renormalization group construction one would try to impose a condition A ~ 1 
by fixing the renormalized coupling Aq of the last renormalization group step to obey that 
condition since in the full construction Aq corresponds to our coupling A. 

The standard nonlinear cr-model has the constraint on the field variable (which we call 
(j) instead of a ) 

02 = 1 . (2) 

Condition (2) can be obtained from (1) by a suitable limit taking K ^^ oo)^ Such a 
constraint however is immediately softened out when starting from the model with a large 
UV cutoff on integrating out high frequency modes, even after the first renormalization group 
step in a renormalization group construction. This can be seen from the renormalization 
group construction of the hierarchical model which has been performed by Gawedzki and 
Kupiainen[4] and later also by Pordt and Reiss[5]. It is rather obvious anyhow: Once you 
have (at least) two independent frequency modes, fiuctuations of one may compensate those 
of the other such that the constraint (2) is restored for the sum. These fiuctuations are not 
even highly improbable since neighbouring frequency modes may look similar in position 
space for frequencies close to the border line between the two. Thus we obtain for K a value 
of order 1 after the first step. The much more difficult part of the ultraviolet analysis of the 
model - so far only performed in the hierarchical case for N > 2 and as long as the effective 
coupling stays small - is to show that the Lagrangian (1) is a good approximation to the full 
model. That implies in particular that the model has only one marginal direction which is 
well represented by the quartic term in (1). So our starting point is reasonable when giving 
credit to the evidence based on the hierarchical approximation. This hierarchical analysis in 
turn agrees with the seminal papers on the model based on perturbation theory by Brezin, 
Le Guillou and Zinn- Justin [6], and the analysis of Brezin and Zinn- Justin [6] also agrees 
with ours on the IR side in the limit N —>■ oo. Furthermore the generally accepted view 
is confirmed by numerical simulations [7] and, which is of great importance in this respect, 
also by the Bethe Ansatz methods based on the exact S'-matrix [3], which show in particular 
that the model has a mass gap. Nevertheless these results are not fully based on well proven 
assumptions and are rather self-consistent than rigorous. So we note that on the other hand 
that doubts against the general wisdom have been raised by Patrascioiu and Seller [8]. 

We take an UV regularized version of (1) as our starting point. The scale is chosen such 
that the UV cutoff A is situated at A = 1. The situation in constructive field theory is 
often complicated by the fact that the expansions around the situation where the degrees of 

^For the analysis of the model in that limit a lattice regularization is probably most appropriate, see also 
[25, 26] and the comments below. 



freedom are to some extent decoupled start from regularized versions which tend to violate 
symmetries of the model in question. The symmetries on the other hand often greatly 
simplify the perturbative analysis if, as is often the case, an invariant regularization for 
perturbation theory is at hand. Fortunately this time we are on the easy side: Once we have 
introduced an interpolating field, which we now call a, the whole analysis of the model can be 
performed without breaking the 0(A^)-symmetry, in complete agreement with the Mermin- 
Wagner-theorem [9], [10]. When the one-component scalar a- field has been introduced we 
may integrate out the 0-field thus obtaining a new interaction given by (the inverse of) a 
Fredholm determinant. For the UV cutoff model it is well-defined in finite volume. The 
infinite volume limit is taken in the end, once the cluster and Mayer expansions have been 
performed, which allow to divide out the divergent vacuum functional. The analysis of the 
Fredholm determinant proceeds similarly as that of the corresponding determinant in the 
case of the Gross-Neveu-Model [1] . It is simplified in the same way as the expansions are since 
we do not have to distinguish different zones characterized by the mean value of the a-field 
- apart from the small field/large field splitting. The main new problem lies in the fact that 
for the inverted Fredholm determinant some of the estimates used to bound the determinant 
(together with antisymmetric tensor products generated by taking derivatives when cluster 
expanding, see [l],p.l69 and more generally [11]) are no more valid. The problem is solved 
by deriving new bounds on inverted Fredholm determinants - in the last part of Ch.3, to 
show stability, by introducing a finer splitting of the large field configurations before cluster 
expanding to make sure that the cluster expansion derivatives always produce small terms, 
and by evaluating the expansion derivatives through Cauchy formulae. 

The paper is organized as follows: Our specific choices for the regulators and the basic 
definitions are presented in Ch.2. They are dictated by technical simplicity. In Ch.3 we per- 
form the small/large field splitting and develop the bounds on the various terms in the action 
ensuing from that splitting, as well as on the non-local operator kernels appearing. In par- 
ticular we show that all the kernels appearing fall off exponentially in the small field region. 
In Ch.4 the cluster-expansion is performed which then allows to control the thermodynamic 
limit and to prove the exponential fall-off of the (two-point) correlation function(s). 

After submitting this paper we learned about two important references on the subject. 
First the author was not aware of Kupiainen's work [25]. Secondly, few weeks after submis- 
sion there appeared a preprint by Ito and Tamura [26] . We close the introduction by shortly 
commenting on these papers. Kupiainen regards the N component nonlinear a-model on a 
unit width lattice for arbitrary dimensions d. He shows that the 1/A^-expansion is asymp- 

^This important and beautiful contribution to constructive physics is maybe not as well known as it 
should be to those working in the field. In part this might be due to its title. 



totic above the spherical model critical temperature Ts, which is zero for d = 2. He also 
proves the existence of a mass gap for these temperatures and N sufficiently large. Without 
attaching much importance to the numerical side we just say what W sufficiently large ' 
means. We read from [25] (see equ.(19)) for the two-dimensional case that for given inverse 
temperature /3 one needs 

Since (3 is to be identified with the inverse coupling 1/A in our language this is basically the 
same as our bound: We require 

since the small factor per cluster expansion step (see end of sect. 4. 4) has to beat the factor 
0(m~^) from the spatial integration per link. Similarly the authors of [26] state their result 
in Theorem 24 for 

and (3 large. They regard the same model as Kupiainen, the N component nonlinear a- 
model on a unit width lattice, for d = 2. Thus [25] and [26] analyse the lattice version of 
(1) where the limit i^ — > oo has been taken, i.e. the Heisenberg model. The result [26] 
only concerns the free energy or partition function which is shown to be analytic in (3, given 
A^ as above. Correlation functions have not yet been treated. It seems clear however that 
their method of proof which, as ours, is based on a small/large field cluster expansion is well 
adapted for that case too. We prove exponential fall-off of the two-point function, extension 
to any connected n-point function is straightforward using the Mayer expansion formulae 
for those, see e.g. [19]. The change in sect. 4. 5 would consist in singling out a connecting 
tree now for n external points instead of two. Kupiainen's result on the other hand is based 
on reflection positivity in the form of chess board estimates. It is not clear how the result 
on the exponential fall-off can be extended to general connected functions in this context, 
so strictly speaking (as he does) his result only holds for those correlation functions which 
have no nontrivial truncations. ^ An important point shared by [25] and [26] (in fact the 
authors of [26] could have referred themselves to [25] here) is that they both apply the 
Brydges-Federbush random walk representation to show and use exponential fall-off of the 
lattice kernels of l/[p^ + rn^ +icr]. In the continuum we only succeed in proving exponential 
fall-off for small flelds a. This is the main reason why we introduce a whole hierarchy of 
large fleld regions with larger and larger protection corridors (see (60)- (63)), and the fall-off 
over the corridors has to make up for the (possibly) absent fall-off in the large field domain. 
Apart from this [26] is technically closer to my paper than to [25]. It is more detailed on 

•^In special cases he succeeds in performing truncations by a clever use of certain Ward identities. 



some aspects of the expansions. A number of bounds take similar form here and in [26]. In 
[26] the building blocks of the cluster expansion are taken to be large also in the small field 
region. This has technical advantages, on the other hand treating many degrees of freedom 
as a whole generally tends to deteriorate the numerical bounds. 

2 Presentation and Rigorous Definition of the Regu- 
larized Model 

We want to show that the UV regularized large A^ a-model is massive, i.e. that the cor- 
relation functions decay exponentially. In our explicit representation we will restrict to the 
two-point function, generalizations to arbitrary 2A^-point functions being obvious. Thus 
formally we study the following object: 



Here Dcp indicates the product of (ill-defined) Lebesgue measures D(f)i, . . . ,D(f)N- Before 
giving sense to this expression mathematically by imposing suitable regulators we want to 
introduce the interpolating field a as announced. We rewrite (3) as 



up to a global field-independent normalization factor. Now we can perform the Gaussian 
integrations over the 0-fields to obtain 

S2{x,y) ^ /D(t(^^)(x,i/) det-^/'(p' + 2(T) e-5w/-' + ^/- (5) 



p'^ + ta 



k2 , A'2 



again up to a global field-independent normalization factor and on rescaling 
{N/2\)(l? . (^T^)(a;,y) denotes the position space kernel of the operator ^4^. Its exis- 
tence in £^(IR^) say, will be clear once the cutoffs and thus the support of the measure are 
specified below. ^ As regards notation we will generally use the same letters for position 
and momentum space objects. This lack in precision in our eyes is overcompensated by the 
gain in suggestive shortness. For the same reason and on the basis of the previous remarks 



'^ When studying higher order coorrelation functions it is preferable to work in the space ^^ £^(]R^) 
and to suppress the exponent N of det instead, because in this case the factor replacing {^q:j^){x,y) will 
depend on the flavour indices. We will adopt this convention only in the last part of the paper where it 
somewhat shortens the notation. 



on the size of the constants appearing we will abbreviate by 0(1) sums of products of N- 
independent constants the largest of which appearing will actually be l/m^ (21)- Without 
making this explicit we pay some attention not to collect astronomic numbers into 0(1). 
By performing a translation of the field variable a according to 

t' = a + im^ (6) 

we finally arrive at 

J p^ + m^ + igr p'^ + m^ 

This time the change of normalization stems from three sources: from the translation, from 
a change of normalization of the Fredholm determinant and from a rescaling of the r'-field: 
t' ^ T = \/tt7 t'- In (7) we introduced the coupling constant 



The value of the translation parameter m is fixed below ((17)-(21)) by a gap equation. This 
eliminates the term in the interaction exponential which is linear in r, and this in turn is a 
prerequisite in the 1/A^-expansion, since that term has a coefficient ~ y/N. Before specifying 
the UV and IR regularizations we note that from the point of view of mathematical purity 
it would have been preferable to introduce them from the beginning. This however would 
have blown up the previous manipulations without a real gain since (3) and (7) are in fact 
to be viewed on equal footing as starting points: They both produce the same perturbation 
theory in 1/A^. 

We now introduce the following regularizations: 
UVl: We set the cutoff scale to be 1 and replace 



p'^p'=P'e' (9) 



in (7). 

UV2: We also introduce an UV cutoff for the r-field. When tracing this back to the original 

interaction (1) it amounts to smoothing out the pointlike quartic 0'^-interaction. To the 

1 r 2 
expression /Dre~2 J "^ in (7) corresponds in rigorous notation integration with respect to 

the Gaussian measure dfis{T) with mean zero and covariance Cs{x — y) = S{x — y), or in 

momentum space Cs{p) = l(p)- We replace the (5-function by a regularized version 



1 + f{p) 



where 7r(p) is defined below, see (28). It is a smooth nonnegative function depending on 
p^ only, bounded above by a constant of order l/rn^ (see (29)). f{p) also is a smooth 
nonnegative function depending on p^ only. It vanishes in the origin, grows monotonically 
with p^ such that a{p'^Y < f{p) < A{p^Y with suitable < a < A < oo, and fulfills 

)(a;-y) = 0, ii\x-y\>l. (11) 



^1 + / 

The last condition is the most important one. That all conditions are mutually compatible 
is rather credible. A proof is in the elementary Lemma 1 in [1].^ There a suitable / (which 
in [1] is further restricted by demanding that it should vanish of high order in the origin) is 
constructed explicitly, basically by starting from the characteristic function of the unit ball 
in position space IR^ and taking linear combinations of rescaled convolutions thereof. We 
should note that it is by no means crucial to have a cutoff with these particular properties. 
Only sufficient fall-off of 1/(1 + /) in momentum and position space are required. So 
e.g. 1/(1 + e^ ) would do. The compact support property (11) is however helpful when 
fixing the final covariance of the model, taking into account large field constraints, see 
(67). It eliminates further small correction terms of similar nature as those appearing in 
5C^ (70), cf. the remark after (79). In short the UV cutoff on the r field replaces the 
ultralocal covariance 5{x — y) of this field by a smoothed compact support version of the 5- 
function sandwiched between the two a/1 + vr-factors. The growth properties of f{p) restrict 
the support properties of the corresponding Gaussian measure dnjir) to (real) continuous 
functions [12] and therefore we need not regularize expressions such as r^ etc. In general we 
will view r as an element of the real Hilbert space £^(1R^, IR). 

IR: As an intermediate IR regularization to be taken away in the end we also introduce a 
finite volume cutoff. To be definite we choose a square 

A C IR^ (12) 

centered at the origin with volume 

|A| = 4n^ » 1, n G IN. (13) 

We then restrict the support of the r-field to A. But we do not restrict the Gaussian measure 
to A from the beginning, because this again would increase the number of correction terms 
later when we perform a configuration dependent change of covariance. We want to avoid 



^ In fact we did not prove monotonicity in [1]. This can however be achieved by a shght extension of 
the proof. We do not include it since monotonicity is not needed here, it might however be useful when 
performing a renormalization group construction on the same basis. 



this, but nevertheless want to suppress contributions in the measure supported outside A. 
We therefore introduce a term 

exp(-i? / r^) , R»l (14) 

in the functional integral, and take the limit R -^ oo later on. Note that absorbing this term 
in the measure and taking the limit right away would amount to restricting the covariance 
to A from the beginning [13]. Again our particular choices for the IR cutoff are convenient, 
but not crucial. 

With these preparations we now obtain the following rigorous expression for the regular- 
ized normalized two-point function: 

SU^,y) = i fdfif{T)i I )ix,y) (15) 

X det-^/2(l + —i— -z^rxA)e-^/:a^-A^' e^^^^+v^^^A^ . 

Preg + ^ 

The partition function Z^ is given by 

Z^ = /rf/xKr)det-^/2(l + —i-— ^^rxA)e-^/>H^-A^' e^^^^+v^^-^A^ . (16) 
J pi^g + m^ 

Here xa is the sharp characteristic function of the set A in position space. Instead of xx 
we will mostly use Px to denote the orthogonal projector on the subspace of functions sup- 
ported in X. From the bounds on the action is given in the next section it is clear that Z^ 
will not vanish in finite volume (see (118)). In the following we will mostly suppress explicit 
reference to the regulators by reg and xk for shortness. 

As announced the value of m is fixed by imposing a gap equation eliminating the linear 
term in r from the action, i.e. we demand: 

When evaluating the Tr, the term /^^^r factorizes an both sides of (17), and we obtain the 
relation 

' J (27r)2 p2ep2 j^ ^2 XK 2\' ^ ' 

For a sharp cutoff at p^ = i ^e would find from this 

(19) 



1 


■ln(. 


1+ rn' 


-) = 
8 


^/^ - 4 



with the solution 

m^ = e-- (1 + 0(— -e--)) . (20) 

Ail 

For the case of an exponential cutoff the integral cannot be evaluated analytically, but it is 
easy to find suitable upper and lower bounds saying that 



47r 



m = CraC ^ , (2f) 

where the constant Cm is close to one (lies between 0.9 and 1.1) for A < 1. For definiteness 
we will assume from now on 

2/7r < A < TT so that e~^° < m < 1/6 . (22) 

Taking into account the constraint (17) we thus obtain for the two-point function 

SHx,y) = -^ [df^f{T){ \ )(x,i/)det,-^/^(l + ^^.^rx)e-^/--A-' (23) 

where we used the standard definition 

det„+i(l + K) = det(l + K) e-^rK^TrKH-M-ir^TrK- (24) 

for any traceclass operator K and n G IN. In an expansion based on the parameter 1/A^ the 
canonical choice of covariance is such that it contains all terms of the action quadratic in 
the field r, possibly up to terms which are suppressed for A^ -^ oo. This is not yet the case 
for (23) since the term quadratic in r from det is not suppressed for A^ large: It contains 
1/A^ from g"^ and N from det^ '^ giving A^*^ altogether. Thus the appropriate presentation 
of the two-point function is rather 

S.»(...) = i,/*c(r,(^;,^^l^-)(..,)det-'^l + _^,,.)e-/..-.'' , (25) 
A corresponding change of definition has also been introduced when passing from Z^ to Z^ 



A = /'^n^r-r^rlpfr^/2. 1 .^^^ A,-«L 



^2 



In (25), (26) dfic{T) represents the Gaussian measure with covariance 

C = {1 + f + PattPa)"' . (27) 

Pa is the orthogonal projector onto the subspace C^{A) of £^(1R^), and n is the quadratic 
part in r from det. In momentum space it is given as 

<P) = -IT- / .^ X2 2 ^ 2 ^ ^ ^2 ^ 2 > 0- (28) 

2 J {2,Tx) q + m"^ [p + qY + m^ 



2 

Since the integral is UV convergent, it is largely independent of the cutoff functions e^ which 
we did not write exphcitly. We find in particular 

-(0) = ^^, (29) 

where C-,^ is again a constant close to 1. Furthermore one easily realizes that 

TT < 7r(0) in the operator sense, or in momentum space 7r(0) — 7r(p) > . (30) 

This can either be done by direct calculation or by noting that 

(r,7r(0)r) =^Tr{yV*) > jTrV' = (r,7rr), (31) 

for the operator 

Vir) = -^-gr. (32) 

By (r, vrr) we denote the scalar product, which is given by J t{x) it{x — y) T{y) in position 
space. Note that V has real expectation values in the real Hilbert space >C^(A, IR). 
For later use we collect the following facts about the operator n and (some functions of) the 
kernels of vr and of l/(p^e^ + m^) in position space. 

Lemma 1: 

a) The operator vr fulfills: < vr < 7r(0) . 

b)The kernel of vr in position space denoted as Tc{x — y) (using translation invariance) satisfies: 

i) \7i{x-y)\ < 0(l)e-2H^-J/l, 

ii) furthermore |a/1 + vr (x — y)\ < 0(1) e"^'"'^"^' for x 7^ y. 

c) < [l/(p^eP^ +m^)](x — y) < 0(l)exp{— m|x — yl} . 

Proof: The proof of a) was given previously. The statement b)i) follows from standard 
analyticity arguments: tt{p) is analytic in momentum space for {ImpY ^ 4r7i^ as is directly 
seen from the integrand in (28) by shifting the integration variable q by p/2. The main 
reason to choose the analytic regulator function e^ was that it does preserve (and even 
slightly enlarge)^ this analyticity domain so that b)i) follows. Coming now to the statement 
in b)ii) we first note that the condition x ^ y eliminates the (5-distribution contribution to 
the kernel so that we may regard in fact 



{VT+^r^ - 1 (33) 

As compared to b)i) we now also have to verify that the real part of 1 + vr stays positive for 
{ImpY < Aw? (so as to exclude a cut, i.e. a violation of analyticity due to the square root). 

^by a factor of |0(1)to^|, see also [1] 

10 



Again explicit calculation simply reveals this to be the case, where the regularization again 
slightly improves the situation. 

c) This statement was proven in [1], Lemma 5. The lower bound follows from the represen- 
tations 



e 



-p 



.2 1 V^/ 1 2 1 -r)2 



-Ti^ m^—^e-^r (34) 



p2 gp2 _|_ ^2 p2 + 1 ^Q p^ + 1 P^ + 1 



and 



1 2 1 -p2 

— r/i — r e 



-^(1 - m'e) + m^e C ds 6''^^'+'^ . (35) 

+ 1 JO 



p2 + 1 p^ + 1 p^ 

Since m^e < 1 it becomes now obvious by explicit calculation of the Fourier transforms that 
the kernel of — — ^ — j is pointwise positive. f\ m 



3 Small/Large Field Decomposition and Bounds 

The representation of the correlation functions according to (25), (26) is well-suited for an 
expansion in 1/A^, since the remnants of the action left in deta are all suppressed by factors of 
1/yN or smaller. We then have to bound the contributions from dets for large values of the 
field variable r to show that it is integrable with respect to the Gaussian measure d^ci^)- 
From our starting point we presume that this should be possible, since there the action was 
manifestly integrable. However, to obtain a convergent expansion of the correlation functions 
we have to perform a cluster expansion which makes visible the decoupling of the degrees of 
freedom with increasing separation in space. The cluster expansion interpolation formulae 
modify all nonlocal kernels of the theory, the modification being different for the measure 
and dets. Therefore one global bound is not sufficient. What we rather need are local 
bounds per degree of freedom. The solution we adopt is similar as in [1] , with simplifications 
due to the fact that we only have one phase, and complications due to the fact that the 
model is not fermionic in origin. The latter implies that certain sign cancellations due to 
the Pauli principle are absent in the outcome of the cluster expansion and necessitates finer 
distinctions on the size of the r-field than in [1] . 

We distinguish between small and (a series of) large field configurations depending on 
the size of /^ r^, where A is any (closed) unit square in A with integer valued lower left 
corner coordinates {ni,n2). Then we sum over the possible choices for all squares. For a 
given configuration we take the union of large field squares, enlarge this region by adding all 
squares below some finite distance from those and divide (roughly speaking) the enlarged 

''This fact will be useful later (see in particular (173)), but it is not crucial. 
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region into its connected components. In the interior of any such component we do not 
introduce interpolation parameters, it is even reasonable not to absorb the quadratic part 
of det2 in the covariance there. Rather we use the large field criteria and certain bounds 
on inverted Fredholm determinants to show that these regions are suppressed in probability 
per large field square A and according to the size of /^ r^. Then the expansion is largely 
restricted to the small field region, where the integrability of dets is assured due to the 
small field criterion anyway. As usual such a cluster expansion with constraints goes in hand 
with a certain amount of combinatorics and technicalities coming from all sorts of correction 
terms. These are controlled by means of the large value of N. We are now going to make 
this reasoning precise. 

We subdivide the volume A into the 4n^ unit squares A specified above and regard some 
given r G £^(A). We say that A G A is a large field square w.r.t. r if 

XK [ T^ > A^^/6 , (36) 

J A 

and A G A is a small field square w.r.t. r if 

\K I T^ < N^'^ . (37) 

We introduce a smoothed monotonic step function 9 G C°^(1R) fulfilling 

«(x) = I » ■ f°' ^ 2 -'" (38) 

^ ^ I 1 , for x > 1/4 ^ ' 

Then we also introduce A-n? factors of 1 into the functional integral according to 



^ Aj^IItaII^ - 1) + fi - gf ^-^II^A 



lA = ^(^^^I^ - 1) + (1 - ^(^^^i:^ - 1)) • (39) 

In (39) we set as usual 

l|rA||^=/r^ (40) 



Now the first factor is decomposed further writing 



'^^Mk _ n = ym/ A-^lrAlla _ 1^ _ fl/ A-^lrAl ' 



n=l n=l 

We then may rewrite (25), (26) as a sum of 2^" terms each carrying for any square A a 
factor which is either the first or the second summand in (39). For a square carrying the first 
factor the functional integral is then split up further according to (41). To fix the language 



12 



we say 

Definition: A square carrying the factor 



OUr) := 1 - 9{^^^^^ - 1) (42) 



is called a small field or s-square. A square A carrying the factor 



(43) 



is a called a large field or 1-square. More specifically we call it an /"-square if it carries a 

factor ^n(||rA||2) resulting from the splitting (41). 

An /-square then only contributes to the functional integral if 



XK\\T^\\l> In'/\ (44) 



an /"-square only if 

and an s-square only contributes, if 



^Nin+m > XK\\ta\\1 > -A^"/^ (45) 



XK\\T^\\l<^N'/\ (46) 

so that we will always assume the respective inequality to hold once a square has been specified 
to be I, /" or s, since in this paper we are only bounding contributions to the functional 
integral. As regards notation we will write Pg , Pi , Pi" and Pa for the orthogonal projectors 
onto functions with support in Ag , A; , A^n and A respectively. Here we denote by A^ C A 
resp. A; C A resp. Ap c A the set of small field resp. large field resp. /"— squares in A. 
Note As U A; = A , UneiN ^/" = ^i ■ Before proceeding further with the 1/s decomposition 
we want to show that the small field condition is sufficient to obtain a small upper bound in 
norm on the operator appearing in det: 

Proposition 2: For r G C^i-^) let Ag C A be a collection of unit squares such that for 
A G As we have 

XK\\ta\\1 < 5/4 iV^/^ . (47) 

Then the operator norm of Ag : C^iA) -^ 'C^(A) satisfies: 

\\As\\ < 0(l)iV-^/i2 < jY-2/5 ^ (4g) 

Here A is the operator Pa 2 \ 2 Qt Pa, and Ag is defined to be Pg A P,. 
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Proof: We first regard A^ = PaAPa for A G A^. For cp E >C^(A) and ||v?||2 = 1 we find: 

|(Aav^,Aav^)| < / / \ipix)rix)Fix-y)xAiy)Fiy-z)TizMz)\ (49) 

< g'F'iO) f \^ix)Tix)TiyMy)\ < g' F\Q) f r\x) < F\0)^N''/' . 

Jx,y JxeA 4 

Here F[x — y) is the pointwise positive kernel (see Lemma 1) 

n^~y) = / ,, ^2 2% ^ 2 ' (50) 

which is obviously bounded by its value at 0, which in turn is bounded by 0(1/A), which we 
absorb in 0(1), which we bound by A^~2/5+5/6^ This proves the assertion for a single square 
A. To go from here to the general case one has to exploit the exponential fall-off of the 
kernel F{x — y) (Lemma 1), which deteriorates the bound by a factor of 0(l/'m^), which we 
absorb in 0(1) and bound it again by A^~2/5+5/6^ g^ ^^^ jg^ ^ ^ 'C^(A) with ||v5||2 = 1- 

\{Asip,Asip)\ < J2 \iAPAip,PA'APA"ip)\ (51) 

A,A',A"eAs 

<0(1)/ E expi-™(^-*(^'^')+'^-*(^'^"))} / \rix)xAixMx)Tiy)xA"iyMy)\ ■ 
a,A',A"gAs ^'y 

By performing first the sum over A" and then over A' and using the bound on r, the Schwarz 

inequality and the fact that (/? is normalized, we obtain the bound 

0{1)N-' Ari/12 V ( / v^^ ) ( / r^) V2 < 0(l)Ar-5/6 ^52) 

This ends the proof. ■ 

As announced we want - for given 1/s-regions - enlarge the 1-regions by security belts 
of sufficient width such that the fall-off of the kernels from Lemma 1 will produce a small 
factor if the kernels have to bridge these belts. This procedure generally will merge together 
some of the different connected components of the 1-region. Let A^^, . . . , A[ be the connected 
components of A^. We say there is a connectivity link between A; and Aj, 1 < i,j < r,i ^ j, 
if there exists some Aj G A^ and some Aj G Aj such that there exists A G A with 

dist{Ai, A) + dist{Aj, A) < 2M , (53) 

where we choose for definiteness 

M = — In AT . (54) 

m 

14 



Then we call h, . . . ,ls the maximal subsets of A^ connected by connectivity links and call 
them connectivity components. Obviously s < r. Now we set 

r = r(/) = {A c A| dist{A,Ai) < M} (55) 

and 

r^ = T{li) = {A c A| Af C /i, dist{A,Af) < M} . (56) 

Thus there is a one-to-one relation between the Fj and the k, and the Fj are connected (in 
the standard sense), and we have 



Ti n Fj- = for i^j, and |J F^ = F . (57) 

1 

In set-theoretic relations we always denote by a set of (standard) Lebesgue measure 0. We 
also introduce the sets •ji which (roughly speaking) lie between /j and Fj : 

n 

7i = {A c IR'I Af c /i , rf^st(A, Af) < M/2 } , 7 = LJTi (58) 

1 

so that 

rf2st(7, A - F) > M/2 - V2 . (59) 

Note that for technical reasons we have defined 7 as a subset of H^, not necessarily of A. We 
do so because this definition of 7 is useful when fixing the covariance in the presence of large 
field configurations such that it has good positivity and fall-off properties (see (67)- (70) and 
Lemma 3). 

The previous definitions now are extended to the situation where we split up further the 
A;-region into the components A^n. If the size of the field is very large we also need very 
large security belts to protect our large field regions - such that the decay of the kernel across 
this belt again assures a small contribution. We start again from the connected components 
A;^, . . . , A[ of A; and say that there is an e- connectivity link (or extended connectivity link) 
between AJ and A;, 1 < i,j < r,i 7^ j, if there exists some Aj G A; fl Aj„/ and some 
A,- G Af n A,n" such that there exists A G A with 

dist{Ai,A) + dist{Aj,A) < {n' + n")M . (60) 



*It requires some (elementary) work to really give an explicit proof of that fact, which amounts basically 
to transferring the square A constituting the connectivity link between A^ and A^- to the centre of a line 
of minimal length connecting A.^ and Aj and showing that then either this transferred square or two of its 
neighbours touching each other connect together A^ and A^ within some Tk- We skip the proof since it is 
not crucial for us that the F, are connected. 
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The e-connectivity components are then the maximal subsets of A^ connected by e-connecti- 
vity hnks. We call them /^, 1 < « < s', and obviously s' < s < r . Now we set 

r = r(/) = U{A C A| dist{A,Ain) <nM} (61) 

n 

and 

r, = r{l',) = U{ A C A| Af„ := A/n n Af C It , dtst{A, Af„) <nM} . (62) 

n 

Again there is a one-to-one relation between the F^ and the If, and as before 

s' 

Tf n q = for i^j, and (J T," = r . (63) 

1 

Starting from the 1/s- decomposition of the volume A we now decompose the Fredholm 
determinant, define the s-dependent final covariance and bound the large field action. With 
the definition of the operator A (Proposition 2) we can write the Fredholm determinant as 
det(l -|- iA). We first separate As from the rest of A via the standard relation 

det-Vl + iA) = det-Vl + iAs) det-^l + ^iA") (64) 

1 + iAs 

with 

A" := A- As = A' + Ai, A' := PsAPi + PiAPs. (65) 

Since A has real spectrum, the operator l/{l + iA) is well-defined. For Ag we now proceed as 
indicated before (see (25)), i.e. we absorb the quadratic part in r into the covariance. When 
doing so we obtain the following (transitory) expression for the inverse propagator C;~^: 

Cf,' = PsTiPs + 1 + /. (66) 

We express (66) in terms of C^^ (67), the basic reason for this being the fact that we are 
not able to deduce suitable fall-off properties in position space for the inverse of (66). Our 
final choice for the configuration dependent covariance will rather be 

C~^ = VTT^ii - P^ + eP^ + f)VTT^ . (67) 

Here e is introduced so that C^ is bounded also in the large field region. We fix it as 

e = N-^o . (68) 

Choosing (67) we have to control the difference between (66) and (67), since it is (66) which 
is isolated from the action. Writing 

C-' - Cfs' = SC^ - Pi (69) 
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we obtain for (5C^ the sum of terms : 



5C, = ^5^(7), (70) 



5Ci(7) = - Ps{ViT^)P^{ViT^)Ps 



5C2{i) = PKvr+^)(i - p^){ViT^)Ps + Ps(vr+^)(i - p^){Vi + ^)Pi 



+ Pi{ViT^){i - P^){Vi + ^)Pi , 



5C3(7) = (1 - PA)(v/rT^)(i - p,){ViT^)p^ + PA(v/rT^)(i - p-,)(v/rT^)(i - Pa) 

+ (1 - PA)(v/rT^)(l - P,){VlT^){l - Pa) , 

5Ci{-i) = ViT^sP^VTT^ . 

Having introduced the final covariance we may now rewrite the expression for the two-point 
function based on the Gaussian measure d/i^ with covariance C^ normahzed such that 



dfi^ir) = 1. (71) 

Since our covariance is configuration dependent there will be a change of normalization of 
the functional integral when changing the 1/s- assignment. Relative to the situation where 
7 = this normalization factor is given by [13] 

Z^ = det^/2(C^/Co) , (72) 

where Cq is given below (76). Taking into account this factor we may rewrite (25) as 

J2i,s Zy J dn^ir) ( o :j2:i„^^J ix,y)G^ 
S^ix, y) = — ^7;^"^.^?:^ , 7 = 7(/) . (73) 



For the action G^ we find collecting the results of the previous manipulations 

l + iA 



G, = Hasa, ^i(r) Hasa. e^ir) e'^'^^-i "' detg ^^'(1 + zA,) det2-^/'(l + -r^iJ^') (74) 



Note that we get indeed dets (1 + iA^^ det2 (1 + ^^^j^ iA") after using the gap equation 
and absorbing the quadratic part of det~ '^(1 + iAg) since 

TrA = TrA, + TrAi, TrAi = Tr{ \ Ai) , (75) 

1 ~r ?'-/ie 
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on using TrA' = and Ag Ai = 0. 

We first analyse the covariance C^. Then we bound the normahzation factors Z^ and 
the correction terms 6C^. Finally we bound the large field determinant. Calling Co the 
covariance C^ for the case that 7 = which means 

we may write the inverse of (67) as 

oo 

c, = Co + Co{c^' - c:;')c, = Co Y.[(Co' - c:;')CoY (77) 

r=0 
1 1 °° 1 1 11 



VTT^ VTT7 ^oW^^TJ '' ' VTTJ' VTTJ VT+ 
1 1 ^ „ , . . 1 „ , . ,„ 1 1 



TT 



^°+vTTTTT7?/'''-^'lTT7^^<^-^'ri 



VTTHi + ft; " "1 + / " " i + fVi 



TX 



The sums are obviously norm-convergent. At this stage the support properties of 1/(1 + /) 
(11) become very helpful. They imply that in position space C^ may be written in terms of 
a simple sum over disconnected pieces with support restricted to (a neighbourhood of) 7j. 
We obtain 

n 

C^ = Co + C^, C^ := ^C^^ , (78) 

i=\ 

If we had only imposed exponential fall-off for 1/(1 + /) , arbitrarily many terms coupling 
the various 7, would have appeared. They could be shown to be small using the distance 
between the various 7^ of size ~ M and the fall-off of 1/(1 + /) and of } , but still they 
would be a nuisance. The fall-off properties of Cq have been analysed in Lemma 1. The 
complications stemming from nonempty 7 are controlled easily in 
Lemma 3: The kernel C"' for 7 7^ satisfies the following estimates: 

\C'^{x,y)\ < 0(l)iV2/5exp{-2m(rfist(x,7) + ciist(i/,7))} . (80) 

For x, y G A — r or a; G Fj , y eTj with i ^ j we find: 

1 , . 



\C\x,y)\ < 0(1) 



JYI8/5 



and for xGF, i/GA — F 



\C\x,y)\ <0(1)^. (82) 
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Finally we have 

|Co(a;,t/)| < 0(l)exp{-2m|x-y|} . (83) 

Proof: We have to control the contribution of the infinite sum over r in (79). We abbreviate 
O = P,. {l-e) £[7^ ^7. (1 -^)r , B = -^ -^ (84) 

r.=0 ■■■ + / Vl + TT i + / 

SO that 

CT. = BOB* (85) 

Obviously ||0|| < N'^/^ and ||i?|| < 1. Furthermore the kernel of B is continuous and 
pointwise bounded by 0(1). By inserting characteristic functions of squares A between B 
and O and between O and B*, summing over the squares, using the fall-off properties of the 
kernels and bounding 

\{XA, OxA')l < N'^' (86) 

we then arrive at the bounds stated in Lemma 3. For the required properties of the kernels 
see Lemma 1 and (11). The minimal distances of points fulfilling the conditions specified in 
Lemma 3 follow from the definitions (54)- (59). ■ 

We remark that the bounds in Lemma 3 could be somewhat improved on by using methods 
similar to those employed in the proof of Lemma 4. We do not do so because this improve- 
ment would not strengthen our final bounds anyway. Note in particular that the cluster 
expansion will be performed such that only C'''-terms bridging the gap between 7 and A — F 
will be produced. Now we are going to bound the factors Z^. 

Lemma 4: Let I7I denote the volume of 7. Then 

1 < Z^ < e°(^)l^l . (87) 

Proof : Using (76), (78), (79) we have 
Z, = det'/\C,/Co) = detV2(i+ 1 ^C^«) = deti/^(l + Y. [(1-^)^7.7^1") (88) 

= n, det^/2(i + E [(1 - ^) Pf. T^D = n. det'/'icjCo) = n, z,^ . 

Again we used the support properties of 1/(1 -|- /) to factorize the determinant. For Z^^ we 
now find 

Z,. = det^/^(l + ^|^~'_^^"'^^_^ ) = det-^/^(l - (1 - .) P,^ ^) (89) 
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= expTr{(-l/2)ln(l-(l-5)P,,^)} = expTr{l/2X: ^(1 - ^^ [P,, ^]^'} . 

This expression implies Z^ > 1. On the other hand we may use Lemma 3' from [1], which 
says that for an Hermitian trace class operator A and orthogonal Projector P we have the 
inequality : 

Tr{PAPy < TrPA^'P . (90) 

Applying this to -^ and Py. we may bound 

T^^J^f P^^y = Tr{Pr. yIj P^J (91) 

<Tr{P,^{^rP,;) <0(1)|7.|. 
using the fact thatQ 

Using this we obtain 

Z, < exp{0(l)|7| E^T^a-^n < exp(0(l)|7|) • (93) 

r>l ^ 

This proves Lemma 3. h 

We now come to the bounds on the correction terms (5Cj(7) from (70). 
Lemma 5: 

i) (5Ci(7) < (as an operator), 

ii) <^C3(7) < 1 + 7r(0), 5^4(7) < 0(l)A^-2/5 (as operators), 

iii) ||5C2(7)|| < 0(l)A^-2, \5C2{l){x - y)\ < 0(1) inf{e-2"^l^-^l, N~^} . 

Proof : i) is immediately obvious from the positivity of tt. 
ii) The first statement is obvious since 

503(7) = iVTT^)ii - Py)iVTT^) - PAiVTT^)ii - Py)iViT^)PA . (94) 

Note that 503(7) o^ily enters through interactions with field configurations of support out- 
side A, which will be suppressed anyway when taking R -^ 00, (Prop. 8, (114)). The bound 
on 504(7) follows from the definition of e in (68). 

^Unfortunately the factor of r~^/^ appearing in (92) is falsely written as 2^*" in [1]. This mistake fortu- 
nately is of no consequence however. 
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iii) The first statement in iii) follows from the exponential fall-off of \J\ + vr (Lemma 1) and 
the fact that dist{Ai, (JR^ — l)) > ^^ (see (58)). This implies a bound on 5C2{'j) of the 
form in iii), closer inspection shows that 0(1) is basically given by m^^, two powers coming 
from the integration over the kernel bridging the distance gap and one coming from a norm 
bound on the second Vl + vr . The second statement in iii) also follows from the definition 



of 7 and from the fall-off of the kernel of ^1 + n . m 

Now we come to the bound on the nondiagonal term det2 (1 + j^ ,^^ iA") in the action 
(74). We need to get a suitable bound for this term which is sufficiently stable under the 
modifications caused by the cluster expansion parameters. We (temporarily) introduce the 
operator B through 

B = , \ , tA" = —-5- — t(^ + A*s)A" ■ (95) 

1 + iAs {l + tAs)il-iA*)^ "^ ^ ^ 

Using the facts that the A-operators have real expectation values in real Hilbert space, that 
TrA^A" = and cyclicity we find 

|det^^(l + 5)| = |det"^(l + 5)| = |det"^(H- 5*)| = det"^/=^(l + D) , (96) 

where 

D = B + B* + B*B . (97) 

Now we may apply the norm bound on As from Proposition 2 to realize that B coincides 
with iA" up to small corrections, more precisely: 
Lemma 6: For ip G jC^{A.) we find 

B(f = iA"(f + 6A"<p, B*<p = -iA"*(f + A"*6*<p , (98) 

where the operator S is bounded in norm as 

< {1 + a)\\As\\ < {1 + a) N~^/^ « 1 (99) 



with suitable < a « 1. 
Proof: Since we have 

the statements of the Lemma follow directly from Proposition 2, where a may be chosen to 
obey an upper bound of size ~ ||As||. ■ 
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Now we can also bound the operator D. For ip G iC?{^) normalized to one we obtain 

(V9, D^) = 2(v9, A'V) - 2(A'V,v9) + (v?, SA"^) + (5A'V, V?) (101) 



Since the first two terms drop out, this entails 

- n\\A"Ln\\n > - 

251-7] 



{V,D^) > {l-v)\\A"^\\l-v\\A"^\\, > -4t^ (102) 



with the choice 

7] = 2\\5\\ «l . (103) 

Splitting the selfadjoint operator D into its negative part D^ and its nonnegative part D^ 

D = D+ - D_ (104) 

we thus have obtained that 

< Z^„ < — ^^ < iV-^/5 . (105) 

- - 251- r] - ^ ' 

Using this we may now proceed to a bound on |det2^"'^(l + B)\ = det""'^'^(l + D). We find 

det-i/2(^ ^ ^) ^ ^-i/2TrD det-i/2(i + D) (106) 

= e-i/^^'-^*^ det2''/'(l + D+) det2"^/'(l-D_) < det2"^/'(l - D_) . 

Here we used again the fact that TrB = iTrAi = —TrB*. Evaluating the trace of D_ in 
an eigenbasis of D_ one may easily establish the bound 

det2-'/'(l - D_) < exp(l/2Tr( ^_^^^^-_^ )) < exp(^-^^^TrD!) . (107) 

To verify the first inequality one observes that for x > — £ > — 1 we have x — ln(l+x) < ^J^ . 
Now it remains to bound TrD^ . We call {<f-} a suitable set of normalized eigenfunctions 
oi D_ and find 

TrD'_ = E(V^-, D'V-) = E[(v^-> D^-)? (108) 

< E(g')ll-'''V-ll)' < ^^'fTr(A"-A") < ^l,fTr(A-A) 



^ i"'*U^''(/i^(s-^'') ^ om-"'9Hif . / 



_2\ 



(2vr)2>2^3 + m2' 
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In the first inequality in (108) we made use of (102). It is admittedly pedantic to insist on 

factors as (4/5)^ in our context. To pass from Tr{A"*A") to the expression in the last line 

in (108) it is sufficient to take away the projectors Pi or Pg and thus to bound Tr{A"*A") 

in terms of Tr{A*A) which is given by the double integral. We have obtained 

Lemma 7: 

|det2-^/^(l + -^^A")\ < e«(i)^-^^^ L^' . (109) 

^ ^ 1 + iAs 

Finally we may also bound somewhat further the terms (r, 6Cs{'y)r) and (r, (5C4(7)r) 

appearing in G^. Writing f := rxm.'^-A we have 

1/2 (r, 5^3(7) r) = (f, 5^3(7) r) + 1/2 (f, 5Csi^)f) . (110) 

We may then bound 

(f, 5^3(7) r) + 1/2 (f, 5Csi^)r) <R/2 (f, f) + ^(^Xa, txa) (HI) 

for R large enough, using Lemma 5. As for 5(74(7) we find 

(r, 6C,i^)r) < 0(l)iV-2/5[(rxA, txa) + (r, f)] . (112) 

Now we dispose of a complete control of the action G^ from (74) and may collect our in 

findings in 

Proposition 8: For R large enough we have 

Z, G,{r) < e-^4^%°«^-^^^/A.^^ e-^/^/H-A-^ . (113) 

Thus we now take the limit i? — *> 00 and absorb the term e Jir^-a in the covariance, 
which implies that we may replace 

C^ ^ PkG^Pa (114) 

and restrict the action to configurations t{x) with suppr C A . 

Proof: The proof concerning (114) is to be found e.g. in [13], so we have only to gather 

the pieces for the proof of (113). In the first term e ^°° ^i we collected together the 

-1/2 f T^ 

contribution e ^i from (74) the A^-contribution from (109), the term from the bound 
on Z^ in Lemma 4, where we used 

0(1) |7| < 0(1) i^^r |A,| < 0(1) (^)2 Ar-V6 f ^2 (115) 

m m JA; 

and the contributions from (111) and (112) in A;. Finally we also absorbed in this term 
a contribution coming from 502(7). Iii the term e -'as ^ we have absorbed the 
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contribution in A^ from (109), (111), (112) and again a contribution coming from SC2{'y). 
Finally we absorbed the one from the bound on |det3^ (1 + ^^s)| to be derived now: 
|det3 (1 + iAs)\ can be bounded using the inequality 

|TrA"| < \\A''-^\\Tr{A*A) (116) 

valid for any traceclass operator A and n > 2. To bound ||A"^^|| we use Proposition 2. 
We may restrict to n = 3, the subsequent terms in the expansion of dets (1 + iAg) being 
much smaller 

\TrAl\ < \\As\\TriA:As) < iV-^/SiV-i 0(l)(rx., 7r(0) rx.) < Oil)N~'/' f r' . (117) 

With the help of the previous remarks and this relation we can verify the bound 

Idetg ' {l + tAs)\ < e ^ ^ -'As . ■ 

The following result now is immediate. 

Corollary: Reducing the volume A to a single square A equipped with a small field condition 
(42) we find 

Z^ = / ^^^(r)^^ = 1 + o(iV-^/^) , (118) 



where G^ is the integrand from (26) restricted to the single small field square volume, and 

dji^ij) is the normalized measure with covariance XaC'oXa- 

The statement follows from the bound (113) restricted to one small field square in A. 

It will be useful later on to bound the large field contribution in (113), r.h.s. by a product 
of suppression factors in probability per square A G F. If (43) holds we may write 



/ r^ > 1/2 / T^ + ^ivV6 and / r^ > 1/2 / r^ + ^Ar"/^ . (119) 



'AeA; JAeA; 8 JAeAjn JAeA; 

Therefore we obtain 
Lemma 9: 

e-^4^^ < e~^/^4^^ e-^^^^l^^l U e"^"^ . (120) 

Proof: It suffices to observe that for A^ sufficiently large we have (using (42), (45)) 

0(i)(!^M)-2Ar«/6 > Art . 
m 

m 
Now we have sufficient control of the action to start with the expansions. 
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4 The Expansions, Proof of Mass Generation 

4.1 The General form of the Expansions 

The cluster expansion allows to control the spatial correlations of the model. When combined 
with a subsequent Mayer expansion, which frees the clusters from their hard core constraints, 
it allows to take the thermodynamic limit and to bound the decay of the correlation functions. 
We proceed similarly as in [1] and use the general formalism for cluster expansions presented 
in [18], which in turn is an elaboration on a theme which has been the subject of several 
seminal papers by Brydges and collaborators over more than a decade. We apply in particular 
the Brydges-Kennedy formulae [15]. For general references on cluster expansions see also 
[13], where the presentation is close to the original way of introducing cluster expansions in 
constructive field theory, and [14], [19], which are close to our way of presentation. 

The cluster expansion is a technique to select explicit connections between different spa- 
tial regions. The best formulas for the clusters involve trees, which are the minimal way 
to connect abstract objects together. We call the subsequent formulae forest formulae, the 
forests generally consisting of several disconnected trees. The basic building blocks of our 
expansion are the large field blocks F^ (62) composed of (generally many) large field squares 
and their security belts, and the individual small field squares A from 

S:=A-r. (121) 

From the point of view of the presentation it seems advantageous to connect together these 
large field blocks F| already by a first cluster expansion, and then to proceed to a second 
one, the building blocks of which are given by the outcome of the first. Then the expansion 
really connects together unit size squares which allows to somewhat unify the language as 
regards convergence criteria etc. In view of the existence of the excellent presentations to 
be found in [14]-[19] and since we stick very closely to [18] we hardly give indications on the 
proofs of cluster expansion formulae here. The general forest formula we are going to use 
will be given now. We introduce the following notation: 

Let J be a finite index set (in our context the set A of the squares A G A) and P{I) the set 
of all unordered pairs {i,j) & I x I, i ^ j. A (unordered) forest JF on / is a subset of P{I) 
which does not contain loops {ii, 12) ... {in,ii)- Any such forest splits as a single union of 
disjoint trees, and it gives also a decomposition of / into |/| — |JF| clusters (some of them 
possibly singletons). The non-trivial clusters are connected by the (non-empty) trees of the 
forest. Let H he a function of variables Xij, ij G P. Then the following forest formula due 
to Brydges is proven in [18]: 
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H{1,...A) = EdI tdh,)({l[^)H){h{^{h)) , (122) 



where 

h^^{h) = ml{huleL^{i,j)} (123) 

and Ljr{i^j) is the unique path in the forest JF connecting i to j. If no such path exists, by 
convention hf^ih) = 0. 

This interpolation formula will subsequently be applied to our expression for the two- 
point function, more precisely to the summands in the numerator and denominator of (73) 
with given //s-assignments. As mentioned we proceed in two steps. The first rather trivial 
one is to connect together the squares in the components F^ of F^. Let Pr<= be the set of all 
pairs (i, j) of distinct squares in F*^. We define 6^ = 0, if Aj fl A^ = or if Aj and A^ belong 
to different components F^(-^^ 7^ ^lij) of F*^, eij = 1 otherwise, and rjij = 1 — Sij. Our first 
forest formula is simply 

1 = E n(^ J' dhi) n iVi + eihf^ih)) . (124) 

This follows directly from the application the forest formula to H{{xij}) = HijePre (^«i^ii + 
rjij) using that here H{1, ..., 1) = 1. 

The only non-zero terms in this formula are those for which the clusters associated to the 
forest J-'i are exactly the set of connected components F^ of the large field region. Indeed 
they cannot be larger because of the factor HieJ^i ^'' ^'^^ '^^^ they be smaller because of 
the factor Yii^j^iivi + ^i^i ^(h)), which is zero if there are some neighbours belonging to the 
same component (for which r]ij = 0) belonging to different clusters (for which hij^{h) = 0). 
Therefore this formula simply associates connecting trees of "neighbour links" to each such 
connected component, but in a symmetric way without arbitrary choices. We remark finally 
that in (124) the interpolated factors Yli^j^j {rji+Sihf (h)) after giving the necessary constraints 
on the clusters can be bounded simply by 1. 

The second cluster expansion links together the previous clusters by interpolating all the 
non-local kernels in the theory. It gives a forest formula which is an extension of the first 
one. We consider all non-local kernels in our theory, that is 



^=^ , VlT^ , ^— ^ . (125) 

Note that due to the support property (11) of 1/(1 -|- /) and our choice of treating each F^ 
as one connected block of the second expansion, we need not interpolate C^ = Cq + C'^ as 
a whole: When all kernels appearing in (125) are interpolated such that they do not connect 
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any more different clusters of the second expansion, then C^ - with these interpolated kernels 
replacing the noninterpolated ones in the expression for C^ - does not connect different 
clusters either. The three kernels from (125) will be generically called K. Now the second 
expansion takes into account the connections built by the first, i.e. it interpolates only the 
links 

Ki{x,y) = K,,{x,y) = A,{x) K{x,y) Aj{y) (126) 

for squares which belong to different clusters of the first forest. Let Z{K, F^, A) be a generic 
name for the quantities we want to compute, namely the numerator and denominator in 
(73). Then the second forest formula gives: 

ZiK, ^^ A) = E n i^i f dhi) n {Vi + eihf{h)) x (127) 

X E n itdh) n {-^)z{K{{h^,.^,}),A), 

where Z{K{\hjr^^j:^\h) is a functional integral with interpolated kernels K{\hj:.^^j:^\ 
These interpolated kernels are defined by -ft' ({/ijc-2_jc-^} ) = hi^' '^{h)Ki{x,y), where hi ^' ^(/i) 
is the inf of the h parameters of the lines of JF2 — JF^ on the unique path in JF2 joining Aj to 
Aj (if / = {i,j)). Again if no such path exists, by convention hi ^' '^(h) = 0. In other words 
the path is computed with the full forest, but only the parameters of the forest JF2 — JF^ are 
taken into account for the interpolated non-local kernels. 

The product Y[ieJ^2-J^i('£r^ ^^ ^ short notation for an operator which derives with respect 
to a parameter xi multiplying Ki where K is any of the non-local kernels, and then takes xi to 
1. Therefore the action of n/GJ^2-j^i(d^) creates the product n/gjP2-.^i ^i (with summation 
over the finite set of possible i^'s), multiplied either by functional derivatives hooked to 
both ends (for the case where the derivatives apply to the measure and are evaluated by 
partial integration) or by other terms descended from action exponential, if the derivatives 
apply directly to the action. In section IV.4 we give the list of the corresponding derived 
"vertices" produced by these derivatives. The important fact to be shown is that because 
these derivatives act on terms which a carry a factor N~^ , a; > , in fact to each such 
vertex, hence to each link of this second expansion, is associated a factor which tends to zero 
as A^ — » 00. 

It is an important property of the forest formulas of this type that they preserve positivity 
properties [18], so that if fi' is a positive operator, K{{hjr^-jr^}) is also positive. This is not 
obvious at first sight form the infimum rule of (123), but it is true because for any ordering 
of the h parameters (say hi < ... < hn) there is a way {which varies with the ordering) to 
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rewrite the interpolated K{h) as an explicit sum of positive operators [18]: 

K{h) = Y.{K - Vi) E Xv,,Kxv,, (128) 

The functions Xv,q ^^^ ^^^ characteristic functions of the clusters built with the part of the 
forest made of lines j9, p + l,...,n. For us (as for anyone interpolating Gaussian measures) 
this preservation of positivity is crucial when the covariance C^ is interpolated. 

4.2 The Cluster amplitudes. Factorization 

From (127) we realize that the quantities Z{K^ F^, A) factorize over contributions, the mu- 
tually disjoint supports of which - to be called polymers - are the blocks connected together 
by the links of the disjoint trees in the forest J^2- So they take the form 

A{K, F^ Y)= J2 n i^i f dhi) n iVi + eihj' (h)) x (129) 

trees {ri4=:ri 'eTi -^o i^T^ 

X E n ifdhi) n {^)A{K{{hr,^rA),Y) , 

trees Ta on Y,T2DTi ^eTa-Ti -^O ieTa-Ti "-^i 

The trees Tia join together the connected subsets of F fl F^, their union, called 7^, (which 
in fact is a forest) becomes a subset of a single tree when adding the links from 7^ — 7^. 
The trees T2 connect together all of the polymer Y, so they have |y| — 1 elements. Then 
(similarly as above (122)) A{K {{hr2-Ti}) , ^) is a functional integral with interpolated kernels 
K{{h'r2-Tx])- These kernels are defined by K{{h'r2-Tx]) = ^«^' ^{h)Ki{x,y), where /i;^' ^(/i) 
is the inf of the h parameters of the lines oi T2 — Ti on the unique path in T2 joining Aj to 
Aj for / = {i,j). 

Now regarding more explicitly the two-point function (73) we get the following formula 
as result of the cluster expansion: 

EiUaZ,^ E A^iYux,y)il/iq~iy.)Ut2A'iYi) 



g,y 



qAf X _ y,ny,=o,u,y,=A 



g,y.' 



with the following explanations : 

1) The amplitudes for the polymers depend on the choice / of the large field region. By 
shorthand notation / stands for the infinite series of possible choices s, l^, P, . . . . Corre- 
spondingly the sum J2i stands for the infinite sum over those choices. We note already that 
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there is no convergence problem associated with this infinite sum due to the suppression 
factors (120). 

2) The difference between the numerator and the denominator in (130) is that in the nu- 
merator there is one external polymer depending on the source points x and y. Note that 
there is no nonzero contribution in which the points x and y lie in two distinct polymers. 
This would necessitate to cut the factor ( i.jn'i+i r — )(''''^) ^^^'-^ ^ product of two pieces of 
disjoint supportj^, one containing x and the other y. Such a contribution obviously vanishes. 
The absence of such a contribution can be traced back to the symmetry — > —0 of the 
action (1). 

Since by the rule of our cluster expansion, each component 7^ of the large field region is 
contained in exactly one polymer Y , we may absorb each normalization factor Z^^ into its 
cluster, defining 

A{Y):=A{Y) n Z,^- (131) 

ahaCY 

The simplest cluster is a single small field square A C 5* = A — V^^ Due to (118) we find 
in this case 

v4o(A) = l + o(iV-i/5) . (132) 

Therefore it is convenient to cancel out the background of trivial single square small field 
clusters, hence to introduce for a polymer Y the normalized amplitude 

Then we obtain the usual dilute polymer representation: 

T.I E a\Y,,x,y)j^Y{Ua\Y,) 

t 
Y^r\Yj^O 

*^<"' "^ = E, E (i/sOnL.oWO ■ ''''* 

To get factorization we must analyze how the choice of / affects the cluster amplitudes. 
The choice of the large field regions A^n for fixed n is a local one, which means that the 
constraints implied by the choice are of finite range. The sum over these choices therefore 



^We write 

1 1 1 



p2 + m2 + igr 1 + ^^^igr p^ + m? 

and interpolate the kernel l/(p^ + m?), see also the proof of the Theorem below . 

^^We assume the square not to contain the external points x, y which may be thought to lie far apart. 
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can be absorbed into the value of (redefined) factorized amplitudes. Indeed we can replace 
the global sums over s, l^, /^, . . . by local ones: 



1 ^jr-^ - 1 ' 



J2 E a\Y^,x,y)J——l[a\Y^= ^ KY^x, y)-——l[biY^ , (135) 

I q,yi W ^)- i=2 q,Yi \H -"-J- j=2 



1 



E E (i/90n«'(>^0= E iVq^mm) (136) 

with the explanations: 

(i) The right sum is over all sets {Yi, ..., V^} where the Yi are sets of A's, a single A being 

excluded, (except if it is an external square containing one of the source points x and y). 

One has the disjointness or hard core constraints Fj fl Y^- = for i ^ j. 

(ii) h{Y) is computed from a{Y) through 

h{Y) = E' a\Y) , (137) 

where the sum is over all assignments of large field regions included in Y. This sum J2' is 
submitted to constraints (as indicated): We define Ai(Y) := AjflF = \JnAin(Y) , Ain(Y) := 
AinClY and sum over the s, /"-assignments within Y with the following restriction: For given 
Y any assignment for which there exists some A G Ap (Y) with 

dist {A,{dY-dA)) < nM (138) 

is forbidden, because otherwise our polymer would not contain the whole of the large field 
block r^ containing A and associated with Ai{Y). It is also evident that it does contain 
this block if (138) does not hold for any square from Ain[Y). With this definition of the 
amplitudes b{Y) we now obtain factorization: 

E KY^, X, y)j-^nU2biY^ 



Q.Yi 



E (l/g!)nLi^(^i) 



4.3 The Mayer Expansion and the Convergence Criterion 

(139) has now the form required for the application of the Mayer expansion in a standard 
way. The hard core interaction between two clusters or polymers X, F is ViX^ F) = 
if X n y = 0, and K(X, F) = +00 if X fl F 7^ 0, and the disjointness constraint for the 

30 



polymers can be replaced by the inclusion of an interaction e ^(^»'^j) between each pair of 
polymers. A configuration M is an ordered sequence of polymers. We define b'^{M) by 

6^(M)=T(M)(in6(r,)), (140) 

where the connectivity factor T{M) is defined using connected Graphs G on M, by 

T{M) := E n (e-^^^"^'^ - 1) • (141) 

G connected on M ij&G 

Then we can divide by the vacuum functional to obtain 

S2{x, y)= J2 b^W , (142) 

M (x, y)— configuration 

where M is a sequence of overlapping polymers Yi, ..., Yg, the first of which contains the 

squares containing x and y and thus includes the factor ( 2 <rr)?+i r — )(-''' ^) fro^ ('''3). The 

sufficient condition for the convergence of (142) in the thermodynamic limit is well known: 

It is a particular bound on the sum over all clusters, containing a fixed square or point to 

break translation invariance [14,18,19]. We state it as 

Proposition 10: 

I E 6(F)el^l| < 1/2 (143) 

y,oey 

for A^ sufficiently large, uniformly in A, \Y\ being the number of squares in Y. 

The fixed point is chosen to be without restriction. For A^ large enough, (143) in fact 
holds if one replaces the number e in (143) by any other constant. To deduce convergence of 
(142) under condition (143) requires to reorganize the connectivity factor T{M) according 
to a tree formula. We can use again the basic forest formula (122) to obtain a symmetric 
sum over all trees. We define 

Vij = (e-^(^-^:') - 1) for i ^ J . (144) 

We call P the set of pairs I < i < j < n. Expanding Y[(ij)ep{^ + %) with (122) we get 
another forest formula, on which we can read the connectivity factor 

TiM) = E IIK.. /' dhi) n (1 + hrit,Mj) , (145) 

where hr{i,j) is the inf of all parameters in the unique path in the tree T joining i to j. 
This formula is then used e.g. like in [14,18,19] to derive the convergence of (142). Remark 
that again every tree coefficient forces the necessary overlaps and is bounded by 1. 
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It remains to prove Proposition 10. We do not give a first principles proof here, but we 
do show how to sufficiently control those contributions to the polymer amplitudes, which do 
not appear in analogous form in e.g. UV-regularized massive y^^-theory, since the latter is 
clearly exposed in many reviews and textbooks, e.g. [14,19,22]. Cluster expansion techniques 
are nowadays applied to much more complicated situations than this, recently also with ac- 
cent on a clear and systematic presentation [20,21]. The aspects not to be encountered in 
a y9'^-treatment are analyzed in section 4.4. Here we reduce the proof to certain bounds on 
functional derivatives generated by the links of the second tree T2 — Ti in (129). Because 
the amplitude h{Y) is given by a tree formula we will sum over all squares in Y by following 
the natural ordering of the tree, from the leaves towards the root, i.e. the particular square 
containing 0. The factorial of the Cayley theorem counting the number of (unordered) trees 
is compensated in the usual way by the symmetry factor 1/|1^|! that one naturally gets 
when summing over all positions of labeled squares [14,19]. Then the only requirements to 
complete the proof of (143) are 

(i) summable decay of the factor associated to each tree link. This is obvious for the Sij links 
of 7i, because these extend only over neighbours, so have bounded range. For the tree links 
of ^2 — 7i, it follows from the decay of the corresponding kernels (125), see Lemmas 1,3. 
(ii) A small factor for each tree link, or equivalently for each square of Y . This will compen- 
sate in particular for the combinatorial factors to choose which term of the action to act on 
by the derivatives etc. For tree links of Ti this small factor comes from the one associated 
to each of the large field squares, hence from Lemma 9. Once a square is chosen large field 
we still have the choices /^, . . . /"°, the value of no depending on the distance of the square 
from the boundary of Y . The sum over the n- values converges (rapidly) due to (120). For 
the tree links of 7^ — 7^ the small factor comes from the negative powers of A^ generated at 
the ends of these links ("vertices"). These small factors are described in more detail in the 
next section. Remark that all types of small factors tend to zero as A^ ^ 00. 
We note that the small factor per square should be there on taking into account the bound 
on the action as a net effect. (113) was derived before performing the cluster expansion. 
Does it still hold once the interpolation parameters and support restrictions are introduced? 
It does indeed, because support restrictions do not cause any harm in the reasoning of Ch.3, 
because all interpolated kernels are bounded in modulus by the modulus of their noninter- 
polated versions (see (128)), and because the interpolated versions of the operator A still 
have real spectrum. Then one easily realizes that all statements go through as before, in 
particular the proof of Proposition 2 and of Lemma 7. A slightly more serious modiffcation 
of the action is caused by the use of the Cauchy formula below, it will be controlled by 
Lemmas 12 and 13. 
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4.4 The Outcome of the Derivatives 

With the tools previously developed we now want to show the existence of the correlation 
functions in the thermodynamic limit. We have at our disposal exponentially decaying 
kernels, a suitable stability bound on the action (Proposition 8), and we have arranged 
things such that derivatives will produce a small factor corresponding to the small coupling. 
As compared to a treatment of UV-regularized 99^. The main new features to be analysed 
are the following: 

a) The action is nonlocal, and the covariance is interpolated twice. 

b) There is a small/large field split, and thus small factors per derivative appear in various 
different forms. 

c) The action is nonpolynomial, which implies in particular that terms descended from the 
action by derivation may be rederived arbitrarily often. 

The amplitudes of the polymers Y are given as sums over trees (129) which are the 
factorized contributions coming from the forest formula (122). When performing the h- 
derivatives those may either apply to dfi^(Y) or to 

){x,y) det^'^/'(l + lAs) det^'^/'(l + —^tA") e^/^C^, 5c,r) _ ^^^q^ 



Prea + ^"^ + 'i'9'TXY 1 + ^^, 



reg 



Here we went back to (73) (remembering that the term e Jd^-a'^ jg now absent, cf. Propo- 
sition 8). In (146) the kernels from (125), which appear in C^ and the action, are to be 
replaced by their /i-dependent versions. We write shortly K{h) for i^d/ira-Ti}) and have 
(see (128), (129)...) 

K{h){x, y) = XY{x)h['''^'{h)K{x,y)xY{y) • (147) 

Application of derivatives with respect to dfi^ is evaluated by partial integration ([13], Chap. 
9): 

d,Jdfi};{h,T) ... = ldfx};{h,T) l^^^-^(d,^c,ihmx - y)^-^ . . . (i48) 

In C^ the kernels S = , "*, are interpolated. Thus dhiC-y{h) is of the form 

dh,C^{h) = {dH,S{h))C^S{h) + S{h)C.,dH,S{h) . (149) 

The supports of the derived kernels, i.e. dhiS{h), are by construction restricted to the two 
squares linked by the hi derivation [18], which adds a link to the previous tree. Therefore 
the r functional derivatives are either directly localized in these squares - in the case where 
dhi applies to the first (second) kernel S{h) in C-^, and we consider the ^ derivative on the 
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left (right), or they are only essentially localized - when e.g. dhi applies to the first (second) 
kernel S{h) in C^, and we consider the -^ derivative on the right (left). In the last case this 
means that the ^ functional derivative is linked to its localization square via the second 
(underived) kernel S{h), which is supported over the polymer in question, see (128), (147). 
It has exponential decay, so the links to squares distant from the localization square rapidly 
decrease with distance. Summing over them gives an additional factor ~ 1/m? . Since this 
tolerable deterioration of the bound per derivative is the only effect of essential localization, 
we may forget about this difference from now on. 

The (7^ — TL)-/i-derivatives can apply also to the terms in (107). To roughly keep track of 
the combinatorial factors involved we note that any /i-derivative may apply to any kernel in 
(146) (~ 10 terms). If it applies to the measure there appear two terms with two functional 
derivatives which again may apply to the action (~ 40 terms). Still one should note that the 
effect of these combinatorics is not very important since going through the terms in detail 
(which we shall not do too explicitly) reveals that most of them give much smaller (in N) 
contributions than the dominating ones. This is also true for the sum over the /-assignments: 
Large field contributions, in particular for n > 1, are tiny corrections due to (120). Therefore 
e.g. all the contributions coming from the terms in 5C^ are unimportant: They are only 
present when A; C 7 is not empty. There is one more source of combinatoric increase of the 
number of terms, namely due to the fact that the derivatives may also act on terms produced 
by previous derivatives. For the polynomial part of the action this may only happen a few 
times. But it needs special discussion when regarding the determinants. So we will now go 
through the various contributions and comment how the derivatives act on them. We can 
be short about 



5C^: In all terms we have the kernels \/l + tt, which fall off as exp(— 2m|a; — y\). The 
contributions are listed in (70). When applying an /i-derivative to (5Ci(7) the small factor 
in N comes from dist{'~f,A — F) > In N/m. Due to the fall-off this gives a factor ~ A^~^. 
We may then e.g. write in the bound for the kernel 

)7T) iTT) 

exp{—2m\x — y\) = exp( — —\x — y\)exp{ — —\x — y\) (150) 

and keep the first factor as a kernel with exponential fall-off and bound the second by N~^^* 
using the support restrictions. This is then the small factor per derivative. Note that we 
could also do without extracting this factor from (150), extracting it as a part of (120) 
instead. The same splitting (150) can be applied when the /i-derivatives act on 56*2(7). For 
5(74(7) we may invoke support restrictions to extract N~^^^ as above, additionally we get a 
factor of £ ~ N^"^^^ . The term 5(73(7) does no more contribute due to the limit R -^ 00. The 
same mechanism produces the small factors also, when we apply the functional derivatives 
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5/5t instead of /i-derivatives. Remember the above remarks concerning essential localization. 
By the derivatives we also produce r-fields (essentially) localized in some square A (two fields 
per h-derivative, one per (5/(5r-derivative). If the square A is in A^, we have the choice to 
perform Gaussian integration or to bound the contribution directly using (46) 

I . . . / Ki{z - X) t{x) K2{x - y) ... I < 0{1)N^/^^ I . . . sup \Ki{z - x) K2{x - y)\ . . . | . 

(151) 
This is maybe the simplest way of doing. Note that in this case we still can keep aside a 
factor of A^~3/4 + 2/i2 ^ A^-1/2 pg^ /^-derivative. If the square is in A;n, the bound is achieved 
using (44), (45) and (120). The above-mentioned rederivation of derived terms allows to 
apply (at most) two 6/6t on an h-derived term so that the factor has to be distributed over 
three derivatives leaving in this worst case A^~^/^ per derivative (without invoking large field 
suppresion) . 

Maybe we should also mention shortly the wellknown and well-solved local factorial 
problem. There is the possibility that a large number of r-fields accumulate in a single 
square A, even when regarding only the polynomial part of the action, namely if the tree in 
question has a large coordination number d at that square: There are d links of the type kj^ , 
u = 1, . . .d in the tree, i referring to A. Then bounding the at most 2d r-fields in A C A; 
using (120) (and the Schwarz inequality) gives 



r 



2l'i^-l/4X,r2 . M 



A 



< AUl (152) 



This is not tolerable as a bound when aiming to prove (143), but the solution is in the fact 
that most of the d squares associated to the links lij^ have to be at a large distance from 
A for large d. Extracting a small fraction rj of the kernel decay we can isolate a factor 
associated to d » 1, which is much smaller than ^. Q For a more thorough discussion of 
the point see [14,19] or also [1]. 

Now we regard the Fredholm determinants. As compared to [1] we have to regard an inverted 
determinant. This is related to the fact that we regard a bosonic model, and it means that 
the sign cancellations appearing as a consequence of the Pauli principle which sometimes 
improve the convergence properties are absent. The inverted determinants are raised to the 
power A^/2. For shortness we will change the notationior the rest of this section and suppress 
this power assuming instead the operators Ag , . . . to act in 0^/^ ^"^i-^)- We assume A^ to 
be even, otherwise we still would have to carry around a power 1/2 (without consequence). 
This change entails that we absorb a factor of A^/2 in Tr as well. 

^^It is of order e~^'^ 
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We rewrite the product of the two Fredholm determinants appearing in terms of a single 
one. This is possible, since the interpolation acts equally on all A-operators. We have 

detg 1(1 + lAs) det2-i(l + --^'^') = det^'(l + ^^) e^''^*^»-i/'(^'"=)'+T+^*^"^ . (153) 

Since the Tr of A" multiplied by any power of Ag vanishes, whereas Tr{As + A") = TrA, 
we may rewrite (153) as 

det^\l + iA) expTr{-l/2{iAgf} . (154) 

The cluster derivatives acting on (154) will be evaluated as Cauchy integrals over suitable 
(large) contours. Similar reasoning has been used by lagolnitzer and Magnen [23] in a 
renormalization group analysis of the Edwards model and earlier by Spencer in the analysis 
of the decay of Bethe-Salpeter kernels [24]. To obtain useful bounds using this method 
requires that the derivatives dh^A are always small in norm. At this stage we therefore 
really need the whole cascade of large field splittings from the previous chapter. We have 
Lemma 11: Let / G 7^ — T^ be a link of the cluster expansion joining two squares A, A' 
such that dhiA = Pa' AP^. Then we have 

1 19,,, All < 0{l)N-^^^^ exp{-mdi} , (155) 

if A is a small field square. Here we set di = dist{A, A'). If A is a large field square in A;n, 
we find 

\\dh,A\\ < 0{l)N-^/\[ T^f/^ exp{-mrfj < 0(1) iVT5-5-2" . (155) 

Proof: The result is obtained in the same way as when proving Proposition 2, if A is a 
small field square. If A is in A;™, the distance between the squares is by our expansion rules 
larger than ^"^^^ which assures (156) through the decay of 21^2 (remember in particular 
(44),(45),(62))" ' " . 

For shortness of notation we introduce 

det'\l + Q) := det-\l + iA) (157) 

and first describe how the derivatives act on (157) instead of (154). Namely we write 

dh, ...dh„det-\l + Q) = (158) 

dai ■■■ <9a„ det~i(l + Q + aidh.Q + . . . + otndh^Q) 



ai,...,a„ =0 
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We evaluate (158) by means of a Cauchy formula for the n independent complex variables 
Oj. The idea is to regain the small factor per derivavtive and the distance decay by choosing 
the a-parameters sufficiently large. We note first that det~^(l + Q + ai^h^Q + . . . + andh„Q) 
is analytic in the a-parameters, see Simon [11], as long as 1 + Q + aidhiQ + . . . + andh„Q 
has no eigenvalues. This restricts the maximal size of the \ai\. We choose the size of the 
arparameter corresponding to the link / as follows: 



Ri ■= \ai\ = N^e—"' . (159) 



We now find 

Lemma 12: If the ai are chosen according to (159) then 

||E«^^^.^II < 0{l)N-^'' (160) 

Proof: For the individual entries in the sum the bound follows on inspection. If the supports 
of the links (i.e. the pairs A, A') are mutually disjoint it stays true, since then the dhiA 
are mutually orthogonal. If they are not, we again employ the argument (see above (152)) 
that in this case the links corresponding to a large coordination number d in the tree have 
to grow longer and longer. In this case the sum may be performed using the remnant decay 

m. J 

e 10 '^'. ■ 

Remark: When proving the exponential decay of the two-point function in the end of the 
paper we would like have exponential decay with mass m up to corrections small with A^ 
(without invoking the analyticity improvement due to the UV cutoff). It may then be 
necessary to use the full decay for at most two links among those appearing at a branch 
point of the respective tree (see below, proof of Theorem). Obviously this does not change 
the norm bound (160) at all, since we may bound the sum in the same way keeping aside a 
fraction of the decay for rf — 2 links only. 

So we now evaluate (158) through 

1^;,, ...a,„det-i(l + Q)| = |(-^)" / 2 ^ ^ det-\l + Q + Y.aidh,Q)\ (161) 



iRi...R„ ai 
^ (t^)" o ^ o snp\det-\l + Q + J2aid,^Q) 

ZTT Ki, . . . , Kn a . 



^^ if these links are indispensable to join via the tree the squares containing the points x and y in the 
external polymer A{Y, x, y). 
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where the sup is to be taken over the arparameters on the circles i?/. Thus we obtain indeed 

per derivative a factor 

— iV-V6e-^'^< . (162) 

27r 

Before ending the discussion of how to evaluate derivatives acting on det we mention how 
we treat the (5/(5rm-derivatives. In this case we choose (in modification of (159)) 

RUr) = N'" . (163) 

We thus collect a smaller factor in A^ from the r-derivative because 5/5t annihilates a 
possibly large r-factor, on the other hand we do not get a distance decay factor and need 
not do so, because it is already present in the term dhiS{h) which accompanies 5/5t (see 
(149)). 

Of course it remains to give suitable bounds on the Fredholm determinants modified by 
the a-parameters. We have to remember that our true object of interest is not det^ (1 + Q) 
but rather the subtracted determinant (154). First we note that we may still evaluate the 
derivatives acting on (154) by introducing a-parameters, on replacing as before for a given 
choice of of h- and r-derivatives 

Aa:= A —^ A + ^aidhiA + ^a^S^^A and similarly for As, A' . (164) 

I m 

So after bounding the Cauchy integrals we have to bound 

sup|det^^(l + iA„) exp{-Tr(A^,)2}|. (165) 

a 2 

The task is to reproduce the bounds on the action from Ch.3 on replacing A -^ A^ ■ In- 
spection shows that the proofs of Propositions 2, Lemmas 6,7 and part of Proposition 8 (as 
far as (117) is concerned) have to be redone with this modification on A. We collect our 
findings in 

Lemma 13: We assume that the kernels A are restricted to a given polymer F C A of the 
cluster expansion. Then we have 

a) \\Aas\\ < 0(l)A^-i/^ (replacing (48)) 

b) \det-\l + j^^JA'^)\ <exp{0(l)iV-i/4/y(r2 + 1)} (replacing (109)) 

c) |Tr(AL)| < V-^/^ IvnA^r' (replacing (117)) 

Remark: Note again that due to our change of notation a factor of A^/2 has been absorbed 
in Tr together with a corresponding change in det. 

Proof: The proof of a) is trivial from Proposition 2 and Lemma 12. As for b) we have 
to go again through the considerations leading from (95) to (109). Since the reasoning is 
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analogous, we will be rather short. Introducing the quantities Ba, Da as we did for Aa we 
find that (95) to (97) stay true also for complex a. The essential modification occurs in 
(101), (102): since the a^ are complex we find 

t{<f, A» - t{^, A'» = ~2j2lmai{<f, dh,A"<f) - 2j2lma^{if, 6r^A"^) (166) 

Oil "m 

instead of for a = 0. Correspondingly we have to modify (102). The norm bound (160) 
then still implies 

Da < 0{1)N~^/^ . (167) 

which is weaker than (105) but sufficient for us. In evaluating TrD'^_ we take into account 
the additional contribution too. Since now 

0< Da- < {2 + r])\6A"\ + ^^2 Ri\dH,A"\ + Y.2 R^r) \6r^A"\ , t] « 1 (168) 

I m 

it is straightforward to realize that we may bound TrD'^_ by 

TrDl_ < 0{l)N'^/\ f T^ + 1) . (169) 

JAnY 

The first contribution is obtained similarly as in (108), it is quadratic in r, but we can keep 
aside a small factor. The additional contribution is proportional to the number of squares 
touched by (5T--derivatives (< |F| ), thus it is independent of the size of r. This ends the 
proof of b). 

c) The proof is as in (116), (117). ■ 

From Lemma 13 we now find that (113) (on restriction to F C A and on using interpolated 
kernels) is to be replaced by 

Z^Glair) < e--^^.-^%°«^"^''/A-.(-^ + i) . (170) 

So the large field suppresion stays unaltered and in the bound on the polymer amplitudes 
there is at most a factor of ~ 1 + 0{N~^^^'^) per small field square from the action to beat 
(we could tolerate 0(1)). 

Here we may end our discussion on the outcome of the derivatives. We have shown that 
we have a small factor ~ N~^^^ per derivative and factor of e~^ per large field square. 
All links are through kernels decaying exponentially with mass > m. This is sufficient to 
beat the factors 0(1) per square from the combinatoric choices and from the action. We 
pointed out that this is sufficient for the proof of Proposition 10. 

^^where we mentioned already that just taking the maximal value gives a crude bound since most terms 
are much smaller than the leading ones 
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4.5 Exponential Decay of the Correlation Functions 

Now we have proven the existence of S2{x,y) in the TD hmit we want to proceed to the 

announced result on its exponential decay. 

Theorem: For N » 1 sufficiently large the inifinite volume two-point function decays 

exponentially 

\S2{x,y)\ < 0(l)e-™'l^~^l (171) 

with 

m' = m(l + o(A^i/i°)) . (172) 

Remarks: 0(1) is an A^-independent positive number. The estimate on the exponent of N 
in (172) is of course not optimal. The proof goes through without much change also for any 
2n-point function. Using the effects of the UV-cutoff we could replace m' by m. 
Proof: The reasoning is very similar to that of [1] though somewhat simpler. The point is 
now to realize that the convergence proof still works when we put aside the decay factor 
appearing in (171). We may assume x and y far apart. They both have to be contained in 
the same polymer A{Y,x,y), and we have to extract the decay factor when calculating its 
amplitude. More specifically we shall extract it from the sum over trees T2 in (129), where 
we ffist only deal with those trees T for which 7i is empty, namely we ffist assume that 
Y does not contain large field squares, which is the dominant contribution. Obviously the 
decay is associated with the factor 

1 11 



p^ + m? + igr ' 1 + -^i^ igr p^ + m? 

which appears in the external polymer. The kernel 21^2 is interpolated and thus in par- 
ticular restricted in support to Y. Let Ai and A2 be the squares in Y containing x and 
y. For given tree T there is a unique path in T connecting Ai and A2. We call it T', 
noting that T' is a tree with coordination numbers di = 2, apart from the ends, where they 
equal 1. Its complement in T will be called T" . It has several connected components in 
general. Each of these connected components may be viewed as being rooted at some square 
(attached to links) from T' . Keeping these squares fixed for the moment and summing over 
the positions of the other squares in the various connected components of T" then provides 
us for these connected components with the usual polymer bound (Proposition 10) sufficient 
for convergence. It remains to sum over the positions of the squares in T' apart from Ai 
and A2, which are sitting on the ends. For given positions of those squares we may isolate 
a factor of 

el^'l n K,{xi,,y,) . (174) 



l'£T' 



40 



Here e ~ o{N^^/^'^) is part of the small factor per small field derivative, the other being 
used to beat the combinatoric constants etc., see above. The kernels Kii{xi>,yi/) are those 
generated by the derivatives of the expansion. They all fall off exponentially with at least 
the rate of ^:^^, so they all may be bounded by the modulus of [-T^;^]{xi',yi') up to a 
constant ~ 0(1), which we absorb in e. The coordinates {xi',yi>) are situated in the two 
squares linked by /' G T' and are to be integrated over those squares.]^ In Lemma 1 we 
showed that the kernel of 21^2 is pointwise positive. From this we then obtain easily that 
(174), when integrated over the intermediate squares and summed over their positions is 
bounded by 

e^'''^[{^^,f\x,y) (175) 

(up to a constant ~ 0(1), which we absorb in e.) Note that having split up the tree T 
does not change the way in which the sum over the trees is performed. We succeeded in 
extracting the factor (175) due to the fact that two squares in the external polymer are fixed 
instead of only one as in Proposition 10. When summing over all possible values of \Y\ and 
using the polymer bound (143) we now obtain a bound of the form 

\S,{x,y)\ < 0(1) [^^(1 + E^'^'' (^^)'^'')](^,?/) • (176) 

pZ _|_ ^Z ^ pZ _^ ^z 

Here the first term is the contribution for |y| = 2 and where the single h- derivative applies 
to the second factor in (173). This is the only case where it does not produce a factor < e. 
Performing the geometric series in (176) now proves (171) on using 

(^- n ){x,y) < 0(1) exp{-{m - e/m)\x - y\} . (177) 

Finally we have to make sure that large field contributions do not spoil our estimate. For 
this it suffices to note that in the large field region we have at our disposal a factor of 
< exp(— A^^/^) per square of P^, half of which may be put aside per each square of P^, on 
which ends some /' G T'. Then we only have to note that this factor is much smaller than 
the factor of e which we loose instead, and that the links within P^ are of short range. ■ 
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